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Abstract. A generalization of the Hopf-Cole transformation and its relation to the Bur- 
gers equation of integer order and the diffusion equation with quadratic nonlinearity are 
discussed. The explicit form of a particular analytical solution is presented. The existence 
of the travelling wave solution and the interaction of nonlocal perturbation are considered. 
The nonlocal generalizations of the one-dimensional diffusion equation with quadratic non- 
linearity and of the Burgers equation are analyzed. 
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1 Introduction 

The classical Hopf-Cole transformation is applied for the solution of the nonlinear diffusion 
equation. There are two well-known nonlinear generalizations of the diffusion equation: with 
quadratic nonlinearity and the Burgers equation. The first of them has applications in dealing 
with plasma and acoustic phenomena [15]. The Burgers equation was initially proposed by 
H. Bateman while modelling the weak viscous liquid motion [4] and later rediscovered by Burgers 
as a simple nonlinear partial differential equation in studies on turbulence [8]. This equation 
can be viewed as a simplified version of the Navier-Stokes equation and related to the heat 
equation via the Hopf-Cole transformation [12,20]. Presently, the number of applications of the 
Burgers equation is immense (see, for instance, [17,41] and references below). 

In the case when the properties of a system in a certain point of configuration or phase space 
depend not only on the properties of this system at this point, but also on the properties of at 
least one point of the environment, we deal with the nonlocal phenomena. As the examples, 
let us mention the well-known prey-predator system of Volterra with delay in ecology [40], 
the ferromagnetic properties of matter in physics [5], viscoelastic phenomena in mechanics [26]. 
From the mathematical point of view, such phenomena are usually described by the integro- 
differential equations [2]. Over the last few years, more attention has been given to a special 
part of the theory of integro-differential equations, the so-called fractional calculus [32,33,37]. 
This approach is applied not only in the theory of fractals and to the above-mentioned, already 
classical, nonlocal phenomena, but also for the description of electrical, biological and diffusion 
phenomena. The latter topic, as follows from the growing number of publications, receives much 
attention [33]. 

Some time ago, two fractional generalizations of the classical diffusion equation were pro- 
posed. One of them leads to replacing the second space derivative by the fractional one 
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where (/> = t), is a fractional derivative in the sense of Caputo [33,37], where < p < 1, 
and a is a parameter of nonlocahty. 

From the physical point of view, we may consider this spatial fractional derivative as a Fourier 
transformation of the fractional power of the wave number k (see a short review in Appendix B). 

From this approach, but for the wave equation, the fractional derivative was considered by 
A.N. Gerasimov in [16]. The other generalization, proposed by R.R. Nigmatullin in [30], is 
related to the fractional substitution of the time derivative 

aDt(t> - a4>xx = 0. 

In this situation, due to the well-known relation between the diffusion and the Burgers 
equations, we may expect two analogous nonlocal generalizations of the Burgers equation 

(Pt + 4^(t>x-aaDl(t) = Q, (1.2) 

and the other one 

aDP<p + ^,Px-act)xx = ^. (1.3) 

These generalizations (1.2) and (1.3) are just the applications of the ideas of [16,30]. Usually, 
equation (1.3) is called the fractional Burgers equation. However, most interesting and perhaps 
most productive is the third, fractional Burgers equation with nonlocal nonlinearity (FBENN), 
i.e. the nonlinear and nonlocal generalization of the diffusion equation, based on a fractional 
generalization of the Hopf and Cole transformation 

cl>t + \ aDl{aDlr^ct>f - a4>x. = 0, (1.4) 

in which <j){x,t), (j)()(x) G M, — oo < x < +oo; t >0 and the parameter a > 0. In this equation, 
aDxcj) = A^~^ ad^4' is the fractional derivative, where A is the length parameter, and adx4' is the 
Caputo fractional derivative (see Appendix B). Thus, all the terms of the FBENN have integer 
physical dimensions. From the mathematical point of view, in the case of linear equation, the 
transition to the dimensionless form is recommended, and in the case of nonlinear equation such 
transition is necessary (for details, see Appendix A). 

A series of exact analytical solutions of this equation, the asymptotic form of the solutions 
and a fractional generalization of the Reynolds number are presented. Concrete examples cor- 
responding to the simplest behavior of fractal solution are analyzed. 



2 Linearization 

The main problem in solving the FBENN (1.4) is its nonlinearity. If for the linear integro- 
differential equation we may apply some powerful methods [34,36], in the nonlinear case we can 
get only general estimates [24]. Therefore, the application of the old method of linearization 
could be helpful in understanding some properties of the FBENN. 

Let (j) he a known solution of the FBENN (1.4), and (p = (j) + eip is a small perturbation of 
solution (j), where ip = ip{x, t). Then ip{x, t), in the first order of e, obeys the linear equation 

i;t + ^xaDl-Pij-aiJxx = 0. (2.1) 

The solution (j) could be obtained by any method presented in the next sections. Thus, ac- 
cording to the general theorems of the uniqueness and existence of the weak sense solution given 
in equation (6.3) below, the solution of the linearized FBENN (2.1) could be found numerically. 

The asymptotic behavior of the linearized solution tp{x,t) at t — )• ±00 is very important for 
the stability of solution (f){x, t). 
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3 Nonlocal perturbations of the local solutions 

Let us take advantage of the fact that in the case of p = and p = 1 the solutions of the 
respective equations are known. To find how the nonlocality changes of the local solution of the 
FBENN (1.4), we shall perform the expansion of the latter two equations in the neighbourhood 
of + e and 1 — e of the parameter p. 

In the neighbourhood of the point p = + e, the fractional derivative has the expansion 

aDP^ix) = 0(x) + eaiV,0(x) + 0{p^), 
here ai^x<t>{x) is the nonlocal operator 

aN.m = -70(x) - r log (x - di, 

J a 

here 7 is the Euler constant. 

Then the expansion of the FBENN (1.4) in the neighbourhood of the respective points is 
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p = + e, (3.1) 

+ (P(l)^-a(l)^^ = {l-e)[aN^(j)'^ -^aN^4>]^, p = l-e. (3.2) 



Thus, as follows from equations (3.1) and (3.2), in the case of weak nonlocality (small values 
of the parameters p or 1 — p) FBENN (1.4) can be interpreted as the respective classical equa- 
tions perturbed by nonlocal terms of a highly specific form. Both these equations can be solved 
numerically. 



4 The travelling wave solution 

Let us change the reference frame and turn to a new variable ^ = x — ut. Then the FBENN (1.4) 
takes the form 

laDl{aDl-Pct>f = a<p" + u<P'. (4.1) 
The solutions of the corresponding equations for p — )• and p — )• 1 are known 
</,(^) = 2< - 2a In [e'^+^'^^^ _ 1] + ^2, 

<A(0 = ^(l-th[^(^ + C)]). (4.2) 

Suppose that according to the "intermediate" character of the evolution equation (4.1), the 
travelling wave solution (f){S^) of this equation is just a fractional "deformation" of the solution 
(^(0) {(,){<f>^^\(,)) with the integer value of the parameter p = {p = 1) 

where aD^~^ fiO = al^'^fiO means a fractional integral (see Appendix B). Indeed, the substi- 
tution of expressions (4.3) into travelling wave equation (4.1) leads to equations for the corre- 
sponding solutions and 4>^^\(,)- 

As a consequence, from (4.3) we get 4>^^\(,) = d^ip^'^^S,) . The assumption (4.3) allows us to 
get the travelling wave solution of the equation (4.1) by substituting the corresponding solutions 
of this equation for the limit cases p — >• and p — )• 1 into expression (4.3). 
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The travelling wave solution of (4.1) is 

n(c + ^) 



2< - 2a log 



exp 



a 



with the asymptotics (\) = —2uc, ^ 1; 

^ T{2-p) (e-a)P 
with the integration constant C, and 

02 - 



p = 0, 



<p <1, 



h + 



i + exp^&^e; 

with the asymptotics — )• +oo) = — )• — oo) = </>2, 02 > 

5 The fractional generalization of the Hopf— Cole transformation 

In the case of the evolutional Burgers equation (BE), 

4>t + <t>(t>x - 0((l)xx = 0, (5.1) 
well known is the simple nonlinear Hopf-Cole transformation 

Wx{x,t) 



P=l, 



(x, t) = —2a- 



(5.2) 



which relates any solution w{x, t) of the diffusion equation (DE) wt = awxx to the solution (j){x, t 
of the BE. 

In the case of the nonlinear diffusion equation (NDE) with quadratic nonlinearity 
1 



•H + -^(Px - (^(t>xx = 0, 



(5.3) 



there is also a nonlinear transformation, which relates the solution w{x, t) to the solution (j){x, t) 
of the NDE (5.3) 



(x,i) = — 2alogtt;(x,t). 



(5.4) 



We see that the above two transformations are the cases of one and the same transformation 
which is as follows 



(x,i) = — 2a a-D^ log w(x, t). 



(5.5) 



for p = (5.4) and p = 1 (5.2), and the transformation (5.5) itself relates the diffusion equation 
solution w{x,t) to the FBENN (1.4) solution 4>{x,t). 

It is more expedient, however, particularly in applications, to use the fractional generalization 
of the Hopf-Cole transformation 

(t){x,t) = -2aaDl\og[b + w{x,t)], 6 G M. (5.6) 

Indeed, by substituting expression (5.6) into the initial FBENN (1.4), we obtain the equation 

Wt - aWxx 



DP 



b + w 



0, 



which turns into an identity when w{x,t) is the solution of the DE. 
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The point is that by changing the scale and the variables we can always obtain the solution of 
the diffusion equation for || 'w{x, t)\\ <^ 1. In this case, the generalization of the fractal Hopf-Cole 
transformation (5.6) becomes even simpler 



(x,t) = —2aaD^w{x,t). 



(5.7) 



From the formula (5.7) in particular it follows that if the diffusion equation solution asymp- 
totically approaches zero, for instance, for w{x,t) x~'^ for x — )• +oo (q > 0), then the FBENN 
solution also approaches zero: w {x,t)-^x-^P+i^ for X — +00, {q > 0). In the case of the 
exponential asymptotic, more convenient is the expression (5.5). 

Let w{x,t) : x G [0, +oo) Ui G [0, +oo). The fractional derivative is the right Caputo 
derivative in the Weyl sense. In the case of the solution of the diffusion equation 



w{x^ t) = exp 



cx (P't 

2a 4a 



the solution of the FBENN is 

(p{x, t) = —2a log 
4>{x,t) = 



cx c t 
a + exp(-- + --5 



p = 0, 



a + exp 



cx I C^t I 

'2« 4o " 



p=l, 



{x,t) 



-2a +ooDl__ log 



cx c t 
a + exp (-- + - 



<p <1. 



(5.8) 
(5.9) 

(5.10) 



Note that solution (5.10) continuously transforms from solution (5.8) into (5.9) when the 
parameter p runs from p = to p = 1. 

It is the fractional generalization of the Hopf-Cole transformation (5.2) that has been used 
to derive the solution (5.10) which interrelates the diffusion and the FBENN solutions. 

Thus, if r(?') is a fractional generalization of the Hopf-Cole transformation, then the inter- 
relation among the nonlinear diffusion equation (NDE), the Burgers equation (BE), the simple 
diffusion equation (DE) and the nonlinear nonlocal diffusion equation can be graphically shown 
as follows 

DE 

rp(0) rp(l) 



i 



NDE f- 



p->0 



FBENN 



-> BE 



6 The interrelation of the solutions 

The existence of the fractional generalization of the Hopf-Cole transformation may produce 
the impression that the properties of the FBENN solutions can be reduced to the nonlocally 
transform solution of the diffusion equation. Below, we shall show this is not the case. 

It is well-known that any solution w{x,t) of the DE wt = awxx under the simple Hopf- 
Cole transformation (f){x^t) = —2aWx{x,t)/w{x,t) turns into the solution <j){x,t) of the BE 
(pt + (fxpx — ci4'xx = 0. The reverse is not true, as the solution (p{x, t) does not obey the DE: if we 
substitute (p{x, t) into the BE, we obtain the equation which will be more general than the DE 

Wt - aWxx = f{t)w, 
here f{t) is any time function. 
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Thus, there exist the BE solutions that can be expressed through solutions of the simple DE, 
and the proper solutions that are devoid of such representation. 

In a similar way, we can classify also the NDE solutions. Some of them are related to DE 
solutions through the transformation (f){x,t) = —2alogw{x,t). However, there exist also the 
proper solutions that are not related to DE solutions. 

It is important that the interrelation of the solutions in the case of the travelling wave 
equation (4.3) can be generalized and is valid for the NDE, FBENN and BE cases 

Cl){x,t) = aDl(t)^^\x,t), (P{x,t) = aDlr^(l)^^\x,t), 

here, as above, (j)'^^\x,t) {4>^^\x,t)) is the NDE (BE) solution for p = (p = 1). This means 
that the solutions of all these three equations, irrespectively of the fractional Hopf-Cole trans- 
formation, are interrelated. Hence, in particular, it follows that 

ct>^'\x,t)=dxct>^'\x,t). 

This interrelation is, in a sense, more general as is valid for both the solutions related to the 
DE solutions and proper solutions. If we substitute the proper solution (/)(^)(x,t) of the BE, we 
will be able to restore the proper solutions of the NDE. 

Let us consider the initial evolution equation (1.4) together with two limit cases at p — )• 
and p — 7- 1 

4>t + ^(t>l - a4>xx = 0, (6.1) 
4>t + <t><t>x - Oi(l)xx = 0. (6.2) 

For p = we have a diffusion equation with quadratic nonlinearity or just the nonlinear 

diffusion equation (NDE) (6.1), for p = 1 - the Burgers equation (BE) (6.2), for < p < 

1 - the integro-differential equation FBENN (1-4) which we may regard as an "intermediate" 

evolution equation whose solutions turn into the solution of equation (6.1) or (6.2) depending 

on \mi (f)^'^\x^t) or Ymi (j)^'P\x,t). 
p-> p— >i 

The solutions to the Cauchy problem (1.4) have to be understood in some weak sense; there 
are several options reported, e.g., in [6,7]. In the context of the present study, let us just say 
that under the weak solution (1.4) we mean the solution of the integral equation 

<P{x,t) = e°*^-,/.o - \ . ^D%Dl-'P(f)\x,t)dT (6.3) 

^ JO 

motivated by the classical Duhamel formula. 

It is clear now that in the same manner we can also classify the nonlocal solutions of the 
FBENN. Some of them, through the fractional Hopf-Cole transformation, are related to the 
simple DE. However, there are also a number of proper solutions which have no such relation to 
the BE. 

7 Evolution of the initial conditions 

Let the initial conditions for FBENN (1.4) and for the diffusion equation be related by the 
expression 
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This aUows us to express the solution (/>(a;, t) of the FBENN through the initial condition (j)o{x) 
and the general form of the solution of the diffusion equation: 

" + 00 I 



{x,t) = -2aaDP 



log ( 1 + / e--^-^'^^y^°^y^dy 



(7.1) 



The expression of the solution of the FBENN (1.4) in the general form (7.1) allows us to 
analyze the time evolution of the nonlocal solution. Indeed, let w{x,t) be the solution of the 
diffusion equation wt = ctWxx- Then the solution of (7.1) can be expressed in the form (j){x, t) = 
-2aaDPlog{l + w{x,t)). 

As follows from the FBENN (1.4), depending on the values of the parameter p we deal with 
not one but with an infinite number or an hierarchy of integro-differential equations. One 
of the most important properties of the FBENN (1.4) is interrelation between nonlinearity 
and nonlocality: for the fractional value of the parameter p we have the nonlinear-nonlocal 
and for the integer positive p only a nonlinear generalization of the Burgers equation. In this 
hierarchy, due to the substitution (t){x,t) — >■ aDx~^4'{x,t), the low-order equations turn into the 
higher-order ones, but in the inverse direction this transformation is multivalued. 



8 Interaction of nonlinear and nonlocal perturbations 

The relation (7.1) between the solutions of the FBENN and diffusion equation allows to consider 
an interaction of nonlocal and nonlinear perturbations. Two or more perturbations moving with 
a different velocity can overtake each other or flow together into a new intensive perturbation. 
The FBENN also describes the interaction process of two or more moving nonlocal perturbations. 
The principle of superposition is not valid for the nonlinear FBENN, but it is valid for the linear 
diffusion equation. The fractional Hopf-Cole transformation (7.1) interrelates the solutions of 
the nonlocal and nonlinear FBENN and of the linear diffusion equation. Thus, if Wi{x,t) are 
the solutions of the diffusion equation, then 4'{x,t) = —2a aDxi^os'^Wi) are the solutions of 
the FBENN. 

For instance, for two solutions of the diffusion equation in the form 



Wi{x, t) = ai exp 



CiX ^ t 

2a 4a 



1,2, 



we obtain a nonlocal and nonlinear interaction of these perturbations: 

p = 0, 
di, < p < 1, 

p = 1. 



(t>{x,t) 
2a 



log {wi + W2), 

1 d r\og[wi[i,t)+W2{i,t)] 



r(i — p) dx 

ClWl + C2W2 
Wi + W2 



{x-cy 



9 The conservation laws 

In the case of the BE (6.2), for x e E^, Vt > 0, (/)(±oo,t) 
conservation value, or the time invariant value inv, of 

r+00 

/(^) = / (j){x,t) dx = iuv, 



since 



dt 



2^ 



dx 



ac 



2^ 



+00 



0. 



bx{^oo,i) = 0, we have a 
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In the applications, this conservation law is called the "mass" conservation law, because if 
t) can be a one- dimensional density or a gradient of any physical, chemical or biological 
magnitude. Then, the I^^^ corresponds to its conservation. 

In the case of the NDE (4.2), if Vt > 0, 0x(ioo,t) = (pxx{ioo,t) = 0, we again deal with the 
conservation value 



/(O) = 0(+oo,t) 



-oo, t) = inv. 



(9.1) 



since by applying the derivative dx to the evolutionary equation (4.2) followed by integration 
we obtain 



ar(o) 



dt 



; dx 



a(px 



1 



0. 



This conservation value shows that the difference in asymptotic values for any time moment 
remains unchanged. If, for instance, we deal with the evolution of the potentials, the conservation 
value I^^^ (9.1) shows that the difference of potentials for x — t- ±00 does not change. In the 
case of the FBENN (4.3), we again deal with a conservation value if Vt > 0, aDx~^ (p{^oo , t) = 

aDl-P(t>{±OO,t)=0 



because 



/+00 
aDlr^4>{x^ t) dx = inv, 
-00 



(9.2) 



dt 



liaDl-P, 



dx 



a 



+00 



0. 



Even this simple example highlights two important properties of the nonlocal conservation 
law (9.2): it interrelates the conservation values of two different dynamical systems, which can 
be of different mathematical nature (e.g., in our case these values are integral and discrete). 

Note that in the "common" case of the nonlocal BE 

<Pt + cf>(l)x-aaDl-P(l) = 0, 
an analogous conservation integral exists at other asymptotic values 

(j){±00,t) = aDl.-P(l){±00,t) = 0. 

At this point, it is not difficult to characterize the "mass" conservation law of the nonlinear 
nonlocal evolution equation 



aaDlS 



0. 



(9.3) 



The magnitude is the invariant of the evolution equation (9.3): 

"+00 



aDl'P(l){x, t) dx 



mv 



for aDx ^</)(=boo, t) = aDx ^^''''^^(ioo, t) = 0. However, in this case the integrability is sac- 
rificed, and the fractional generalization of the Hopf-Cole transformation does not exist. 
The "energy" of travelling excitation 



2 



+00 



,Dl-P<pfdx, 



(9.4) 
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if Vt > 0, ^(/>(±oo,t) = aOl: ^(p{±oo,t) = 0, as in the case of the BE, is not unchangeable, 
but it is constantly decreasing: 




/+00 r+oo 
{aDl-P^) ■ {aDl-P^) dx = -a / {aDl-P(l)) dx < 0. 
-oo J — oo 

Like in the case of the nonlocal "mass" conservation law, the "energy" K (9.4) links the 
energy of the travelling excitation in the case of the BE (p = 1) and the one-dimensional density 
of energy in the case of the NDE (p = 0). 



10 Symmetries 

The conservation laws are obviously related to the group of the FBENN automorphisms. In 
the case of the usual BE, there is a large Lee symmetry group of point transformations, which 
contains the Galilei, dilaton, and projective transformations and is generated by operators 

D = 2tdt + xdx - (pd^, K = t^dt + txdx - {tu + ^x) d^f,, 
Pi = dx, P2 = dt, B = 2tdx - d^. 

Symmetries in the discrete BE have been studied in [18]. 

An important and in the general case uninvestigated problem is to find the nonlocal and non- 
classical symmetries of the FBENN. Some general aspects were already presented in [1]. Promi- 
sing seems an attempt of computer symmetry analysis, as was done for the nonlinear heat 
equation in [11]. 

Note here one property that allows us to get new solutions of the FBENN. Let v{x, t) be 
a known solution of the FBENN (1.4), and u{x,t) is the solution of the linear equation 

ut+ {aDl-Pv)ux-aUxx = 0. (10.1) 

Then 



(p{x, t) = —2a log u + V 

is a new solution of the FBENN (1.4). At an integer value of p we have the local equation (10.1); 
in particular, for p = 1 we obtain a new solution of the classical Burgers equation. 

Note here that the existence of transformation T expressed by relation (5.6) allows to solve 
the problem of the FBENN (1.4) symmetry group. If Gi is a symmetry group of the diffusion 
equation, then G = TGiT^^ is a symmetry group of the FBENN (1.4). 

Here, we shall note an important feature of the FBENN, related to the FBENN symmetry. 
Let us apply the operator dx to the FBENN and take the operator aDx beyond the brackets 



0. 



On substituting the variables while 

x^y±A{t + B), aDP(P{x,t) ^ aDP(t>{y,t) ± A, 

where A and B are constants, the evolution FBENN will not change its form. This way of 
generating the new solutions is particularly efficient not in the case of travelling excitations 
when (p = (i){x — ut), but in the general case when cp = (p{x, t). 
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One of the ways of constructing exact solutions of some of the nonhnear equations consists 
in finding the corresponding Baklund transformations. The Baklund transformations have been 
found for the majority of the nonhnear equations that are integrated by the inverse scattering 
method [14,31]. 

It is possible to show that the generalization of the Hopf-Cole transformation T^^^ is a sepa- 
rate case of the Baklund transformation 



W. + ^ib + w) aDl-P^P = 0, Wt + ^ [{b + W) aDl-PcP] ^ = 0. 



(10.2) 



On removing (p{x,t) from the system we obtain the diffusion equation wt = olWxx- 

To obtain the FBENN from the system (10.2), we shall perform the following procedure. From 

the diffusion equation for the function w{x,t) and the identity [wt/ib + w)]x = [wx/{b + w)]t 

follows the condition 



b + w 



a 



b + w 



(10.3) 



and from the first equation of the system (10.2) it follows that 



1 



and 



LDx 



+ 



b + w 2a" " b + w 2a ' Aa'^ ' 

Substitution of the expressions (10.4) in the condition (10.3) gives us the FBENN (1.4) 



(10.4) 



0. 



The first equation of the system (10.2) is a fractional generalization of the Hopf-Cole trans- 
formation 

(/)(x,t) = -2aaDl\og[b + w{x,t)]. (10.5) 
On substituting expression (10.5) into the FBENN we obtain the following equation 



np 

a^x 



1 



b + w 



{wt - awxx) 



0, 



:io.6) 



from which it follows that any solution of the diffusion equation can be transformed, through the 
generalization (10.5) of the fractional Hopf-Cole transformation, into an FBENN solution. The 
reverse implication would be incorrect, as from equation (10.6) follows a more general diffusion 
equation. For instance, 

Wt - aWxx = (Ci + fit)){b + w). 

It is because of this above-mentioned ambiguity of the fractional generalization of the Hopf- 
Cole transformation that the FBENN solutions are differentiated into those related to the 
diffusion equation and proper ones. 



11 The asymptotic form of solutions 

A very convenient dimensionless quantity which is used in the nonlinear BE is the Reynolds 
number. This number is just a ratio of the nonlinear and the dissipative terms: Re ~ (jycpx/ct'Pxx- 
In the case when Re <C 1, the influence of the nonlinear term is negligible, but for Re 2> 1 
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this term plays a crucial role and leads to the nonlinear Riemann equation which describes the 
simplest type of the shock waves^. 

For the FBENN (1.4) we may introduce a dimensionless generalization of the Reynolds num- 
ber, which at the p = 1 coincides with the classical Reynolds number for the BE 

Re ~ - - 



where A is the characteristic parameter of length in the model. 

Depending on the value of this number, we obtain two limit cases of the FBENN 

(Ai - a(t)xx = 0, Re < 1, (11.1) 



1 



<Pt + ^ aDPiaDl-Pcpy = 0, Re » 1. (11.2) 

The equation (11.1) is just a diffusion equation, and the equation (11.2) we will call the 
fractional Riemann equation. In the case of Re 3> 1, depending on the value of the order of the 
fractional derivative p, we have three evolution equations 



1 



<j)t + -rx = o, p = 0, 

(^t + laDP{aDl-P^)' = 0, 0<p<l, 

(pt + 4'4>x = o, p = i. 

The "mass" conservation law predetermines the asymptotic form of the FBENN solution. To 
confirm such a result, let us consider some estimates. From the general form of the FBENN 
solution (7.1) it follows that the limit t — )• +oo corresponds to a rather low value of the parame- 
ter a. At a low a, to calculate the values of the corresponding integrals we can apply the saddle 
point approximation. 

The critical point yo can be determined from the equation + a-Cy~^'?5'o(yo) = 0. Then 
the asymptotic expression of the FBENN solution acquires the form 



(x, t) — aDx 



2t 



+ C 



tr{3-p) 



+ Ciix-a)-P. (11.3) 



For X —7- +00 and < p < 1, the solution (j){x,t) — )• (x — a)'^~P/[tT{3 — p)]- Thus, we obtain 
a power-deformed perturbation of the usual solution of the Burgers equation. Note here that 
these estimates are valid not only in the environment of the meaning p = 1, but also for any 
p G M. 

We have to show the region of the validity of solution (11.3). In the limit case, the integral 
in the expression of the mass conservation law diverges. Therefore, for x > xq the solution 
(f>{x) = 0. To determine the value xq, we insert the asymptotic form of solution (11.3) in the 
expression of the mass conservation law (9.2). This means that Xg/2t ~ /. Thus, the maximum 
meaning of the solution 



c(x,t)~ ^ and xo~v2It. 

(2t)i 



'^This equation has some names: Hopf, weak shock waves, Riemann-Hopf. We are foUowing reference [35]. 
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12 The supersymmetric nonlinear nonlocal diffusion 
evolution equation 

We shall show that the FBENN has a supersymmtric generalization. Let the superfield x = 
aDlr^4> + '4> unite two fields of different properties: the "bosonic" field 4>{x,t) and its spinor 
superpartner ^(x,t); is the constant Majorana spinor. The transformations of the fields (/>, -0 
are nonlocal because of the fractional derivatives aDxf{x) 

6r,^P = 7]aDl'P(l), 6r,aDl-P^ = 7]xPx. (12.1) 

However, the commutator of the two transformations (12.1) is a spatial translation 

The supersymmtric equation 

Xt = (xx + IxT^x) (12.2) 



(here P = 9dx + de is a supersymmetric derivative) is a system of two evolutionary equations 

i^t = i^.x + ^{aDl-P(t> ■ V'),, (t>t = + laDP [{aDl-P(t>)^ - V'V'x] (12.3) 

which are invariant in respect of supertransformations (12.1). In the general case, the sys- 
tem (12.3) is a system of two nonlinear nonlocal evolution equations, which becomes local when 



p = 



p=l 



^(Pt = (Pxx + ^{<Pl- V'V'x), 

Ipt = V'xx + ^{4>'^)x, 
^(pt = (pxx + li<P^ - V'V'xOa;- 



The supersymmetric equation (12.2) and the corresponding system of equations (12.3) unite 
two fields of different nature, and only one of them is nonlocal. 

Here, the main point is a general note related to the application of the nonlocal systems. 
Suppose the case when a dynamic system is characterized by two interacting fields, of them one, 
for instance, the "fermionic" field 'ip{x,t), is measured in the course of experiment, whereas the 
other, the "bosonic" field (j){x,t), is assessed only phenomenologically. Actually, such assessment 
in the class of local evolution equations may result in a qualitatively erroneous mathematical 
model of a dynamic system. 



13 Conclusions 

It is important to note that the influence of nonlocality can be arbitrarily great. Therefore, we 
do not describe nonlocality by an additional term in the Burgers equation. 

Let us remind here that the classical Burgers equation belongs to a unique group of the three 
completely integrable second-order PDEs. I suggest that the FBENN also belongs to a unique 
group of the completely integrable nonlocal PDEs of the fractional order. 

The fractional diffusion process is related to the non-Gauss statistics; this results in slow 
diffusion correlators ((Ax)^) oc Df^ with 7 7^ 1, and D is a generalized diffusion coefficient of 
the dimension L'^/T'^. In our case, the FBENN is related to the so-called Levi statistics [39]; at 
the same time, the initial Burgers equation as well as the diffusion equation are related to the 
usual Gauss statistics. 

From our point of view, the FBENN has at least two important advantages: 
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i) the influence of nonlocality is not assumed to be insignificant; 

ii) the relation of the FBENN to the usual diffusion equation allows a lot of analytical solu- 
tions of the FBENN. 

Besides, despite the nonlocality in the proposed nonlinear and nonlocal FBENN, 

• space-localized solutions are possible; 

• nonlocal perturbations in a system described by the FBENN can interact; 

• the Reynolds number is a universal dimensionless parameter for both the local and nonlocal 
Burgers equations; 

• there are nonlocal analogies of both the momentum conservation law and kinetic energy 
dissipation. 

In some fields of physics, we need the vectorial form of the Burgers equation, e.g., in as- 
trophysics to describe the large-scale structure of the Universe [28,29,43]. In such cases, the 
vectorial FBENN can be proposed 

<Pt + l aD%Dl-P4>f - aV(V0) = 0, 

where (f) = (0^, . . . , (f)^) G M", aDx is the fractional generalization of the gradient operator V. 

Note also that for a = from FBENN (1.1) follows the fractional generalization of the 
Riemann equation, which also has numerous applications. 

The proposed FBENN, because of its general character, allows a wide range of applications. 
Actually, we may try to introduce the nonlocal generalization in almost all fields where the BE 
is applied. These are the nonlocal effects in shock wave propagation in acoustics, the effective 
model of the process of nonlinear heat distribution in the environment in the presence of heat 
sources and sinks, the Kardar-Parisi-Zhang (KPZ) equation in the crystal growth phenomena in 
(l-l-l)-dimensions [23], the nonlinear dynamics of moving lines [22], formation of galaxies [29,38], 
the behavior of the magnetic flux line in superconductors [21], and spin glasses [13], as well as 
numerous examples of the application of the usual Burgers equation, presented in the above- 
mentioned monographs [17,41]. 



A The normalized form of equations 

In applications, usually considered are both the non-normalized form of the diffusion equation 
with quadratic nonlinearity 

+ a(/)2 - /30^^ = (A.l) 

and the Burgers equation 

(t)t + a4)(t)x - I3(t)xx = 0, (A.2) 

in which the concrete sense of the coefficients a and (3 depends on the content of a model 
under description. In the dimensional form, j3 is normally related to the coefficient of diffusion, 
whereas in the dimensionless form it relates to the inverse Reynolds number Re~^. 

Equations (A.l) and (A.2) in the non-normalized form, as well as the non-normalized equation 
FBENN (1.4) 

ct^t + aaD%Dl-P<l}f - /3cj)xx = 0, 
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on performing the transformations 

(j){x,t) ^ -4>{x,t), t^^i, 
a p 

are reduced to the normahzed form 

(l)t + + (1)^^ = 0, (j)t+ aDl [aDlr^cf) ^ + (/.^^ = 0, 0t + </)(/.^ + = 0. 

Note that in the case of the Burgers equation (for p = 1) there exists a special transformation 
related to the changed scale of the independent variables without changing the function, 

/3 /3 
a 

which allows also a reduction of the Burgers equation to a dimensionless form (5.1). Thus, now 
we may apply any form of the FBENN or its limit cases the BE or the NDE depending on our 
requirements. 



B The fractional calculus 



The left-side Riemann-Liouville fractional derivative of the order < a < 1 is 



d /■* fir) dr 



r(l - a) dt Ja {t - tY 
The right-side fractional derivative is 



Db-fit) 



r(i-a) 



d 

dt 



f{r)dr 
(t - 0° 



The left-side Riemann-Liouville fractional derivative of the arbitrary order a G M is 



Da+fit) 



1 



d 



n nt 



f{r)dT 



T{n-a)\dt) (t-r)i+"-'^' 
The right-side fractional derivative is 



1 



r(n — a) 



d 
dt 



f{r)dT 
t (r-t)i+— ' 



where n = [a] + 1. 

The corresponding regularized fractional derivative of the function f{x), or its Caputo frac- 
tional derivatives (in honour of M. Caputo [9]), are defined as follows 



1 



d /(r)dr /(a) 



r(i - a 

the right-side derivative being 



drJa {t-rY (t-a) 



r(l-a) 



'd^ 



f{r)dT f{b) 



{r-tY (b-ty 



The Caputo fractional derivative of the arbitrary order a G M is 



1 



r(n — a) 



d 



n ft 



f{r)dT 



it-r 



1 l+a—n 



{t - ay 



(B.l) 
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and the right-side derivative is 



r(n — a) 



dr J Jt (r-i)i+"-" (b - ty 



(B.2) 



If f{x) is absolutely continuous on [a,T] or corresponding [T,b], then the left Caputo frac- 
tional derivative is 

the right-side Caputo fractional derivative being 

where a represents the order of the derivative, n — 1 < a < n. 

The forms (B.3) and (B.4) are often used in physical literature. However, the forms (B.l) 
and (B.2) are more useful, since they may be applied to a wider class of functions. 

The Caputo fractional derivative or integral in the Weyl sense is the corresponding Caputo 
operator for the absolutely continuous function /(x) defined on the whole real axes M. 

Some properties of fractional derivatives and integrals are listed below 

D-^ = i:^f{t) {D^^ = lUit)) , a > 0, 
I- /(t) = -i- r 7^^^^^, t > a, 



^ f' fir)d r 
r{a)Jt (r-t) 

D:+f{t) = I-^m, (B.5) 
no^+i f(r\\ 

^(*) = E rnl l-X ^' + Rn{t), n = [Re a] + 1, 
^ r(l + a + j) 

where Rn{t) = i:t' D^X' f it) ■ 

The derivatives of some functions are: 

Df, sin At = A" sin (\t+^\, -ooD?: cos Xt = A" cos f At + 



where A > 0,a > —1; when a < —1, we have to use the property (B.5); 
_oo A+e^*^'' = A"e^*+'^, Re A > 0. 

The Riesz fractional derivative and integral of order a are defined by the Fourier transfor- 
mation 

Ds/(x) = (ifcr(-F/)(fe)) , i-/(x) = {\knTf)ik)) . 

The Riesz fractional integral could be presented as a convolution 
I^/(x)= / K^ix-OfiOdC, «>0, 
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where Ka{x) is the Riesz kernel 



— 7~"^(a)|3;|"~" In a = n + 2k, 



where k £N and coefficients 7„^(a) are 



2"W2r(§)/r(^) , a^n + 2k, 

(_l)(n-a)/22"-iW2r(f )r (1 + ^) , a = n + 2k. 



Note that the correlation of the Riesz and Caputo fractional derivatives in the Weyl sense 

2cos(a7r/2) ' 

Two special functions of those often used in fractional calculus are as follows: 
the Mittag-Leffler function 



CXD ^ 

z 



_^r(an + /3)' 
the generalized exponential function 

" Z^r(l + a + n)" 

The basic aspects of fractional calculus and its applications could be found in [3, 10, 16, 19, 
25,27,28,30,33,37,42]. 

C Supersymmetry 

Let us consider an important example. Two large infinite classes of groups, Osp(A^|M) and 
SU(A^|M), are used in applications. The orthogonal group O(A^) preserves the invariant XiX^, 
and the group Sp(M) retains the invariant OmCmnf^'^, where the Cmn matrices are real antisym- 
metric matrices and 6i are Grassmann-valued. The orthosymplectic group is now defined as the 
group that preserves the sum 

Osp(7V|M) : Xix' + emCmnO'' = inv . 

Note that the orthosymplectic group contains the product 

0=p(iV|M) = ) (C.l) 

Let us write the generators of Osp(l|4) as Ma = {Pfi, M^^i,, Qa) which have the commutation 
relations 

[Ma,Mb]± = f^AsMc- (C.2) 
In the explicit form, the commutators involving the supersymmetry generator are 
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Let us define superspace as the space created by the pair Xn 9a , where 9a is a Grassmann 
number. Let us define the supersymmetry generator as 

Q« = - ii7'0)ad,. 

The anticommutator between two such generators yields a displacement 

Note that eQ makes the following transformations of the superspace 

x^^ x^- iej^9, 9a^9a + 
We can construct the operator 

Da = ^+i{rd^9)a. 

This operator Da anticommutes with the supersymmetry generator, = 0. This re- 

lation is very important because it allows us to place restrictions on the representations of 
supersymmetry without destroying the symmetry. This permits us to extract the irreducible 
representations from the reducible ones. 

Acknowledgements 

The author would like to express his gratitude to Professors B.A. Dubrovin, M. Pavlov and 
L. Alaniya for the invitation and kind hospitality during the Conference "Geometrical Methods 
in Mathematical Physics" (Moscow State University, December 12-17, 2011). 

References 

[1] Abraham-Shrauner B., Guo A., Hidden and nonlocal symmetries of nonlinear differential equations, in Mo- 
dern Group Analysis: Advanced Analytical and Computational Methods in Mathematical Physics (Acireale, 
1992), Editors N.H. Ibragimov, M. Torrisi, A. Valenti, Kluwer Acad. Publ., Dordrecht, 1993, 1-5. 

[2] Agarwal R.P., O'Regan D. (Editors), Integral and integrodifferential equations. Theory, methods and appli- 
cations. Series in Mathematical Analysts and Applications, Vol. 2, Gordon and Breach Science Publishers, 
Amsterdam, 2000. 

[3] Bardou F., Bouchaud J. -P., Aspect A., Cohen- Tannoudji C, Levy statistics and laser cooling. How rare 
events bring atoms to rest, Cambridge University Press, Cambridge, 2002. 

[4] Bateman H., Some recent researches in the motion of fluids. Monthly Weather Rev. 43 (1915), 163-167. 

[5] Bertotti G., Hysteresis in magnetism: for physicists, materials scientists, and engineers. Academic Press, 
New York, 1998. 

[6] Biler P., Funaki T., Woyczynski W.A., Fractal Burgers equations, J. Differential Equations 148 (1998), 
9-46. 

[7] Biler P., Karch G., Woyczynski W.A., Critical nonlinearity exponent and self-similar asymptotics for Levy 
conservation laws, Ann. Inst. H. Poincare Anal. Non Lmeaire 18 (2001), 613-637. 

[8] Burgers J.M., The nonlinear diffusion equation. Asymptotic solutions and statistical problems, D. Reidel 
Publishing Company, Dordrecht, 1974. 

[9] Caputo M., Linear model of dissipation whose Q is almost frequency independent. H, Geophys. J. R. Astro- 
nom. Soc. 13 (1967), 529-539. 

[10] Carpinteri A., Mainardi F. (Editors), Fractals and fractional calculus in continuum mechanics, CISM Courses 
and Lectures, Vol. 378, Springer- Verlag, Vienna, 1997. 



18 



P. Miskinis 



[11] Clarkson P.A., Mansfield E.L., Symmetries of the nonlinear heat equation, in Modern Group Analy- 
sis: Advanced Analytical and Computational Methods in Mathematical Physics (Acireale, 1992), Editors 
N.H. Ibragimov, M. Torrisi, A. Valenti, Kluwer Acad. Publ., Dordrecht, 1993, 155-171. 

[12] Cole J.D., On a quasi-linear parabolic equation occurring in aerodynamics. Quart. Appl. Math. 9 (1951), 
225-236. 

[13] Fisher D.S., Huse D.A., Directed paths in a random potential, Phys. Rev. B 43 (1991), 10728-10742. 

[14] Fokas A.S., Zakharov V.E. (Editors), Important developments in soliton theory. Springer Series in Nonlinear 
Dynamics, Springer- Verlag, Berlin, 1993. 

[15] Galaktionov V.A., Posashkov S.A., New exact solutions of parabolic equations with quadratic nonlinearities, 
U.S.S.R. Comput. Math, and Math. Phys. 29 (1989), 112-119. 

[16] Gerasimov A.N., A generalization of linear laws of deformation and its application to problems of internal 
friction, Prikl. Mat. Meh. 12 (1948), 251-260. 

[17] Gurbatov S.N., Malakhov A.N., Saichev A.I., Nonlinear random waves and turbulence in nondispersive 
media: waves, rays, particles, Nonlinear Science: Theory and Applications, Manchester University Press, 
Manchester, 1991. 

[18] Hernandez Heredero R., Levi D., Winternitz P., Symmetry preserving discretization of the Burgers equation, 
in SIDE III - Symmetries and Integrability of Difference Equations (Sabaudia, 1998), CRM Proc. Lecture 
Notes, Vol. 25, Amer. Math. Soc, Providence, RI, 2000, 197-208. 

[19] Hilfer R. (Editor), Applications of fractional calculus in physics, World Scientific Publishing Co. Inc., River 
Edge, NJ, 2000. 

[20] Hopf E., The partial differential equation Ut + uux ~ f.iUxx, Comm. Pure Appl. Math. 3 (1950), 201-230. 

[21] Hwa T., Nonequilibrium dynamics of driven line liquids, Phys. Rev. Lett. 69 (1992), 1552-1555, cond- 
mat/9206008. 

[22] Hwa T., Kardar M., Avalanches, hydrodynamics, and discharge events in models of sandpiles, Phys. Rev. A 
45 (1992), 7002-7023. 

[23] Kardar M., Parisi G., Zhang Y.-C, Dynamic scaling of growing interfaces, Phys. Rev. Lett. 56 (1986), 
889-892. 

[24] Kilbas A. A., Generalized fractional differential and difference equations, in Fractional Integral Equations, 
Proceedinds AMAD-2003, Belarusian State University, Minsk, 2003, 5-13. 

[25] Klafter J., Shlesinger M.F., Zumofen G., Beyond Brownian motion, Phys. Today 49 (1996), no. 2, 33-39. 

[26] Lakes R.S., Viscoelastic solids, CRC Mechanical Engineering Series, CRC Press, New York, 1998. 

[27] Miller K.S., Ross B., An introduction to the fractional calculus and fractional differential equations, A Wiley- 
Interscience Publication, John Wiley & Sons Inc., New York, 1993. 

[28] Miskinis P., Nonlinear and nonlocal integrable models, Technika, Vilnius, 2003. 

[29] Miskinis P., The nonlinear heterogeneous diffusion equation and distribution of matter in the early Universe, 
Astrophys. J. Lett. 543 (2000), L95-L98. 

[30] Nigmatullin R.R., To the theoretical explanation of the "universal response" , Phys. Stat. Sol. B 123 (1984), 
739-745. 

[31] Novikov S., Manakov S.V., Pitaevskii L.P., Zakharov V.E., Theory of solitons. The inverse scattering method. 
Contemporary Soviet Mathematics, Plenum, New York, 1984. 

[32] Oldham K.B., Spanier J., The fractional calculus. Theory and applications of differentiation and integration 
to arbitrary order. Mathematics in Science and Engineering, Vol. Ill, Academic Press, New York - London, 
1974. 

[33] Podlubny I., Fractional differential equations. An introduction to fractional derivatives, fractional differential 
equations, to methods of their solution and some of their applications. Mathematics in Science and Engi- 
neering, Vol. 198, Academic Press, San Diego, CA, 1999. 

[34] Polyanin A.D., Manzhirov A.V., Handbook of integral equations, CRC Press, Boca Raton, FL, 1998. 

[35] Riemann B., Uebcr die Fortpflanzung ebener Luftwellen von endlicher Schwingungsweite, Abh. Konigl. Ges. 
Wiss. Gdttmgen 8 (1860), 43-65. 



A Generalization of the Hopf-Cole Transformation 



19 



[36] Sakhnovich L.A., Integral equations with difference kernels on finite intervals, Operator Theory: Advances 
and Applications, Vol. 84, Birkhauser Verlag, Basel, 1996. 

[37] Samko S.G., Kilbas A. A., Marichev O.I., Fractional integrals and derivatives. Theory and applications, 
Gordon and Breach Science Publishers, Yverdon, 1993. 

[38] Shandarin S.F., Zel'dovich Ya.B., The large-scale structure of the universe: turbulence, intermittency, struc- 
tures in a self-gravitating medium. Rev. Modern Phys. 61 (1989), 185-220. 

[39] Shlesinger M.F., Zaslavsky G.M., Frisch U. (Editors), Levy flights and related topics in physics, Lecture 
Notes m Physics, Vol. 450, Springer- Verlag, Berlin, 1995. 

[40] Volterra V., Lecons sur la theorie mathematique de la lutte pour la vie, Gauthier-Villars, Paris, 1931. 

[41] Woyczyhski W.A., Burgers-KPZ turbulence. Gottingen lectures. Lecture Notes in Mathematics, Vol. 1700, 
Springer- Verlag, Berlin, 1998. 

[42] Zaslavsky G.M., Chaos, fractional kinetics, and anomalous transport, Phys. Rep. 371 (2002), 461-580. 

[43] Zel'dovich Ya.B., Gravitational instability: an approximate theory for large density perturbations, As- 
tronom. Astrophys. 5 (1970), 84-95. 



